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ABSTRACT 


The Economic Order Quantity (E0Q) model is extended to 
a deterministic, two-echelon model with an arbitrary number 
of activities on the lower echelon. Two variations of the 
model are developed using minimization of time-average cost 
as an objective. A no-stockouts-allowed case is examined, 
and a method for finding the optimal solution is developed. 
A backorders-allowed model is derived and partially solved 
here in general. A full solution is presented for a re- 
stricted range of lower echelon parameter values. Examples 
of the no-stockouts-allowed model are given and solved. 
The solutions from this model are compared LOuEnOSe Genived 
assuming the activities operate wholly independently. 
Digna canempocential reduction in variable time-average 


cost through the use of this model is demonstrated. 
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I. INTRODUCTION 


A field of primary interest in Operations Research, 
especially as applied to supply support problems, is the 
Study of inventory models. Inventory theory is the appli- 
cation of mathematical optimization techniques to the problem 
of deciding when and in what quantities certain goods should 
be purchased to meet the demand of customers while minimizing 
time-average cost [1]. 

The fundamental inventory model, usually called the 
economic order quantity (EOQ) or Wilson model, presumes that 
all costs and demand are fixed over time, that resupply 
occurs after some deterministic leadtime following placement 
of an order, and that demand will be met without incurring 
shortages of stock, all of this for a single activity hold- 


ing only one good or commodity (see the sketch below). 


Purchase 
every T units of time 








Inventory 
Held 






Constant Demand 


More complex and realistic models have been developed by 
relaxing some or all of these assumptions [1]. Nevertheless, 


the majority of the problems that have been solved concern 





finding Optimal policies for single activities since multi- 
echelon problems are more difficult to formulate and solve. 

For many practical applications inventory models for a 
single activity suffice. But for some large retail firms 
and military supply systems which have centralized procure- 
ment for a system of retail outlets or stock-points, a 
mathematical model which accounts for the multi-activity, 
multi-echelon nature of the inventory system is more appro- 
priate. A classic principle of optimization theory is that 
any allocation of resources for a multi-component system 
based on optimal solutions for each component operating 
independently can be no better than the solution found for 
the system as a whole by subordinating the roles of the 
independent components. It is the author's belief that most 
multi-echelon inventory systems tend to operate sub-optimally, 
because the above-mentioned independent component approach 
Hemusedsto solve for stocking policies. 

Hose paper presents a simple, deterministic, multi- 
echelon model which is optimized on a "system" basis. The 
model is an extension of the EOQ with two echelons, as 


depicted below. 


pie. 
ee 
= 





Purchase Qo 


Inventory 
Held 


Order Ord he 
Qi Order Qm 
bw 
Lrnven Cory Tnventory Inventory 
Held Held 2 Held 
* Constant Constant Constant 
Demand D, Demand D Demand Dm 


~ The upper echelon consists of a single activity which could 
represent a centralized procurement and warehousing facility. 
The lower echelon has m independent retail outlets or stock 
points which replenish from the upper echelon and meet the 
demands of the customers outside the system. A minimum 
time-average cost solution for the system, presuming no 
stockouts are allowed and subject to the system cycle length 
being an integer multiple of the cycle length of each lower 
echelon activity (possibly a different multiplier for each 
activity), 1s derived in Chapter II. In Chapter III, an 
algorithm is developed to find the best choice of integer 
multipliers. The case in which customers' demands may be 
backordered is taken up in Chapter V. Some examples of the 


no-stockouts-allowed model solution are presented in Chapter 





IV with an analysis of the cost benefit available in the 
use of this model over the use of independent EOQ models 
EOrU stiemactavities. 

The work in this paper is basically an extension of 
the two-echelon model with one activity at each echelon 
sketched in [2] and presented in greater depth with the 
finite production rate and backorders cases in [3]. 

Several authors have attacked special cases of the multi- 
echelon problem. A probabilistic model for repairables is 
presented in [4]. A probabilistic two-echelon model is 
solved in [5] using a dynamic programming technique. 

The deterministic models in [2] and [3], while they are 
fine in themselves, solve problems of completely specified 
$size. The stochastic models of [3] and [4] are not multi- 
echelon extensions of the probability models presented in 
Chapter Four of [1]. This paper is a natural two-echelon 
extension of inventory theory's fundamental model with an 
arbitrary number of activities at the lower echelon. It 
is hoped that it also provides a starting point for further 
development of more complex deterministic models and more 
general stochastic models such as extensions of those in 
Chapter Four of f1]. 

The author presumes in his presentation that the reader 
is familiar with the EOQ model and its backorders-allowed 


extension. 





II. THE NO-STOCKOUTS-ALLOWED MODEL 


A. MODEL FORMULATION 

siepose that there 1S an inventory system consisting of 
two echelons of inventory, performing the procurement, the 
holding, and distribution of a single good. The upper echelon 
is a single activity which orders and buys the good from an 
outside source and holds it in inventory for further distri- 
bution within the system. The lower echelon has m independent 
activities which order and resupply from the upeer echelon 
and hold inventory to meet the demands of their customers. 
In setting up the cost equation for the system the following 
are assumed: 


1) All orders are filled immediately (a deterministic 
lead time has no effect on the optimal policy) ; 


Zjeeali=demands are met immediately; 


3) the demand rate is constant and continuous over time 
but may be different for each lower echelon activity; 


4) the goods are purchased only by the upper echelon; 
aadeene DrIce 15 COHStant Over time; 


5) ordering costs are constant over time; 


6) inventory holding costs are products of the annual 
interest rate and the purchase price of the good 
(following [1]); and 


7) the cycle length (time between orders) will be constant 
for each activity; and system cycle length will be a 
positive integer multiple of each lower echelon activ- 
ity's cycle length; the number of complete sub-cycles 
for each lower echelon activity may be different. 


10 





The sketch below is a schematic of the material flow in 
this inventory system. The construction of the mathematical 


model of the system parallels that of the EOQ model. 


Qo 
Activity 
© 
Q Qe 
4 ae 
| ™M 
D,* time D,xtime Dux time 


Throughout this paper whenever a summation sign (2) is 
used, it is assumed that sum is taken for the index variable 
ranging from one to m unless it is specifically stated 
Otherwise. 

Activity 0 is the upper echelon activity. Activities one 
through m are the lower echelon activities. The subscript of 
a variable identifies the parameter or variable with its 
associated activity. The Aj; I.,; and D. aie mee SOc @elny cals ae 
the ordering cost, interest rate, and consumer demand rate 


gia 


fOmeene 1 aeeiyvtty, these, walong with €, the unit price 


of the good, are the system parameters. The Q;> T;, and n; 


ata 





are, respectively, the ordering quantity, the cycle length, 
ele 


and the number of complete cycles over i for the i aet hVlty. 
Q; and T,; are continuous variables while n; 1s a positive 
integer variable. 

Let Z be the total costs of running the inventory system 
fou Onemsyseem-cycle, To: SInce Cachmeyvcle 1S «cCOMmpiiscanor 
n; sub-cycles for the ioe 2GmNVv iby. tOtal OLdelmIng NGOs ts same 
Ay + 2,n,Ai . Total purchasing cost is merely CQ)» for the 
goods are bought only at the upper echelon. Let B;(t) be 


tae quantity Of Inventory on hand at time t. Then the total 


™ Te 
holding costs are ra riCc | B; (t) dt ; PUbbingSthece panes 
tegether ; 


To 
Zz = A.t 7, n. Ai + C$ Qo + ATC J Bede. 


A typical realization of the on-hand inventories within the 


system over one cycle is shown in the following sketch. 


Activity 
Oo 


ASIANA ee 


Activity 
m 


?ANAANAAAAANA #922 
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The assumptions of the model result in: 


Ts =n, TT; fom ail *1 ; 
Q. = D. T; PoOlwdl lea. vad 
(2-1) 
Sol aed, “i 
Te. 
Now consider ie { B: (t) dt for 1st4™, Since 


eine n; sub-cycles are identical and the area under the on- 


ana curve fer one sub-cycle is just that under a triangle, 
aT 
fs 1 Quit 
TC [Brdt = nic fBiwde= nic A | 
° 9S 


FOneeme Uppemmechelow activity Bo {t) can be decomposed 
into the sum of the inventories held for each of the m lower 


echelon activities. If Ee (t) is the inventory held for the 


jth echelon, then By (t) = & Bo (t) where 


. + 
B, ()= {ni-d = Int [¥] | Q; and Tnt[X] is 
Piemomearcstmintcoer Less than or @€qual to X. 


Then 
Ve. nit 
(Bie = [ (n-t- Int fe] ade 


= Q. T. 2, (n.~4-5) = z QM nin-L). 


Summing over 1 gives 


To 
! 
(Bude = LE@Trceb, 
It follows that the total costs per cycle are 


2 Ast Dndi+ CQ + TZhCQi+ LETCQnumeiT. 


PS 





B. MINIMIZATION WITH RESPECT TO Q 
To compute the time-average cost for the system, divide 
the cost per cycle by the cycle length, ie. eile T, . iinas 3 


K= ae + = sere? S + Dec + FETC m-1)9@. (2-2) 


Plenouchwtnts 1S a function in Zm + 1 continuous vari- - 
ables, it can be reduced to a function in one continuous 
variable by using the equations (2-1) that arose from the 


model's assumptions. Since Q. ao T; and We = Ticgel. <= 





a i 
n T; EOGe l,jJe=> L,...,m; Equation (2-2) becomes 
AnD, Aj nD, ; TiC D n,Q, ee n.Q,(n;-1) De 
K= nQ, +h “wad * c2Ddi+ 7 aD, ni i 2 Zon. : (2-3) 
himents tenn, K will be referred to as K(Q, ,n, pe+e,N). 
For a given choice of Ny>-++-oN, K(Q,) will represent equation 
i252): 


Taking the derivative of K(Q,) gives 





3k =) = A.D, _ ., An. D, - DD neo.) yes n,(n.-4) D. 
dQ. Qe n, ae 2D = ID ea 
K 
and sto implies 
dQ, '/ 
z 
‘ 2D? {Ac sIn Ai| 
Q. ia C ny {EPL ED (2-5) 
‘ n: an . 
Now, 





2K Zine D, Dy 
d = n, Qe - 4s, Qs rn Ai > O 


| 


a 
& 
"y 


Since all of the parameters and variables are positive. 
A 
Therefore, Q, gives the minimum of K(Q,). 
Nn 
Substituting Q, given by equation (2-5) into equation 


(2-3) gives 
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V2 


o(n.-4 Di. + T.di 
K(Q,)= [ac [Ase Tndi{{ 7 Boe Pee POE OR (2-6) 
Let Onan... Hm ) be an optimal value 
Gri. fn )» Then 


A * »% ~ 
Theorem 1-1: (Omi) = K CQy, 11, bee im pe 


Puoer: Since igen ann ) < min K (QM), 
—— @,,Maycy Mm 


MCG ny.) Im )< LaGGians, 1.) Now, for each fixed 
Nyy Amy, Ky Gin eeanies)i< KiG@ee--, 1, ) for all Qi. 
Then K CQin-,--aaie ) < K MOM ne, nt ) BOs all Qy- 
Tersy CQians,-..Mm)2 K(Qtn",....n2), Finally, 
deme. ne) < AG, ,n*,-.- Mn) $ K(Qint...Ma-) which 
Hieemeres the result. 
+ be % * 
Comoibary loz: Qe @MON 50m). 


Ce eeOoING REMARKS 
* * * 
Gaven Q, ,n poee on, the remaining components of the 
optimal operating doctrine can be derived. Let starred 


variables represent the optimal values of the variables. 


Then 
* \ a 
i > 2 (2-7) 
net : sa ie alli 5 an 
e D: T;* for all i 5 ez =Op 
1 t —-# 
he = Vl, T, 5 and (2-10) 


ILS 





Qe = DL ne an (211) 


In this chapter, the no-stockouts-allowed model has been 
developed, and the values of the continuous variables which 
minimize the time-average cost function (given any set 
i. - 9 Nin of the integer variables) have been found. 
In the next chapter, a method for finding the optimal values 
of the integer variables is derived which will, by using 
equations (2-7) through (2-11), give the optimal policy for 


operating the inventory system. 
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ee NDI NGSOPTLIMAL 1 5...,N,,.0N THE 
NO-STOCKOUTS-ALLOWED MODEL 
Chapter II developed the no-stockouts-allowed model and 
gave in equation (2-6) the minimum value of K for a choice 
Seei--,lge inethis chapter a method for finding optimal 
M)>+++sNp 1S given and, for the case Aj = A, for all i, an 


optimal search algorithm is stated and proven. 


fee tie GENERAL CASE 

eee cmeewMauipmem (Nl) >...,N,) Cf positive integers: 
Pere (i) be the optamaie kk foreach choice of n as given by 
Bewation (2-6). Let n* be the optimal value of n. 

Minimizing K(n) could involve considering a countably 
Pirie wnunber sOredlternatives for n*. Another potential 
PEoolem in seeking n* 15 that there may exist a sequence 
oemn's, call it fn 5 | "such that the Wag) K(n,)=O, Thus, 
Pvenecnough K(n) 15 positive and, hence, bounded below, it 
may not have a minimum. + 

Fortunately, neither of the above problems arise with 
Kinjeeeemeoremn S-1 and its corollary . given below show that 
Pc minederinite, One-Sided neighborhood G of the point 
eoreerwnosc boundary e€an be computed easily. Hence, n* 
can be found by considering a finite number of candidate 


POmMtES Nn. 


Sis 


: As an example of this peculiarity, consider y(n) = 
for n= 1,2,3,... Clearly, y(n) >0 for all n; yet x Yln=o. 
Hence, there does not exist a finite n* such that y(n*) ¢ y(n) 
for all n. 


ey 





Theorem 3-1: For each ";,14j¢m , there exists an \F 
Pueimunat h(n) > K(l,...,1) whenever nj 2 = regardless of 
the choice of Ns» 1 #j- 
Proof: Let - 

Liny= 3¢ { Kon) - CLD: | 

- f ayen,Ac}{z See: Ee 

Note that meni (laos, 10) if and only if 
mea) > KC,... 4). Let Ve 1 41,2,-.ym} be arbitrary 


but fixed. Let 


Xx = Aet Enki ard 
Won. -)ib. fei D; 
Bt ba 


Then a> O and A. >O . Then 


he D;@;-1) J. Dd: 
ee nA; | {RR + i + B. 


(n,-4) A, TDi + Ay — 


Ln) 





+f B, i ag HTS Red I, bs ye. n: Ail 


N:, - 


Let 


M (n,)= (n,-1L) ArT; ie ee 


Since L(n) - Mn) > O, Ltin)> Mon). Let 
b= Glee...) ‘ and let 


Meet pa 1) +4 


where Int [+] is the greatest integer function. 
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Sal Hels M(n,) tomas erictiyerncreasine= function m7 


No» M(n,) > M(N,) for ns > Ni. But 


me A, LD; 


RGN; ) > (BALT:D: | A, T.D, + A, 13D; 


N 


Seen ed de 
Hence, for all n, 2 N; and regardless of the value of 


Me n., ie), Levy> M(n;) > Lct,.-.,1). 


Ae 
icmerouenewsame conaitions on n, K(n) > K(1,...,2); 
and the proof is complete. 

Theorem 2-1 shows that the neighborhood 
G = fn: den. cN: and tsism| 1s a hyper box in 
the positive integer lattice in m-space with one corner at 
the point (1l,...,1l). The following sketch shows what G 


Poems 11ke 1f nm = (Ny, Nyy). 
The 
Hyper box 


a 





LL 
n 
Ni-L N Ny+1 


a 





Miminitcanrate 1mplrcation ©: theorem 3-1 is the folowing 


eououlanye that n* li1es in the bounded hyper box G. 


Comoe. «FOr K given by equation (2-6) 

Me) Ge 

wetter m* = (1,...,1), or it is not. If 
ieee, ssl), tien n* = G. Assume now that 
ieee, lien K(l,...,1) > K(n*). Now 
cor all n & G, 


: Kemiee= KCle 1) > K(n*) 


which implies that n* is not an element of the complement 
SaGweenemcce mms eG. Inverther case n*¢ G. “This 


Gompliekes the proot. 


The previous theorem and corollary show that in the 
general case n* and, hence, the optimal policy (Qo ie ) 
can be found in a finite number of steps by an exhaustive 
search. The exact number of steps required is equal to the 


Pompe: Of POSItive integer lattice points in G; namely, 


“A 

a (Ni-4) ; For example, suppose that m = 5 and 
c=! 

N. = 11 for i=1,...,5. Then the number of iterations 


reaqumred to find n* is 100,000! 


B. THE CASE Ay = Ay FOR ALL 1 

The balance of this chapter is devoted to developing an 
optimal search algorithm to find n* for the special case of 
having the same order@mng cost at each of the lower echelon 


activities; that is, A. 


oe A, for all i. The method of search 


is demonstrated by the following example in which n = (n,,n,). 
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Ooo scuNnm—ao and No, = 3. The search begims by com- 
puting K(1i,1f) { actually L(1,1)} and saving (1,1) as the 
first candidate for the optimal solution. The two adjacent 
lattice points in G, (1,2) and (2,1) are compared to find 
out which one gives the smaller value for K. Say that 
meee) => K(Z,1). Now K(Z,1) is compared to K(1,1) to find 
Omemmaether (2,1) will become the updated candidate for n*. 
If K(2,1) < K(1,1), then (2,1) becomes the updated candidate. 
Whether the candidate changes or not, the search continues 
from the point (2,1). Now the two points adjacent to (2,1) 
LOmewiteh elther ny, OF M4 increases (that is, (3,1) and 
(2,2)) are compared in the same manner as were (1,2) and 
fet). Sapethat K(2,2) < K(3,1). Then K(2Z,2) is compared 
BUC cdnanddte™ Lon n= as before, and either an update 
occurs or it does not. In either case the search continues 
from the point (2,2). The search terminates when either 
ie Or lens Nemeonc ,eat that point, the coordinates of 
the last candidate give n*. The diagram below shows how 


the above search began and how it may have continued. 


Faint on Patu Candidate n* 
cist) Ct, t) 
(2,1) ue), 
G22) C2525 
(3,2) (2,2) 


(3,3) Stop. n*= Cae) 
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The algorithm, itself, will be stated at the end of 
the chapter. Some theoretical groundwork is required to 
prove that this method of search works. 

In the following development, let g(n) be a positive 
increasing function in each n, defined over the positive 
integer lattice in Euclidean m-space. Let We = (ls s:eeerene les 
Tere No.4 = Cope 1) Bema Cleientuan a Sequcnece Of ies 
which track the progress of the search through the positive 


ieeeceer lattice. Then let 


AQ,= MO. x PAs Ret Ld — Gln a) 5 


44es™ 
and let ny be such that 
= = N 
Ad, 4 CM) Te ee. 
Then the sequence [ a, { traces the coordinate-wise, unit 


ep -sized pam ot steepest ascent in g(n). 


Define h (%) by 


Thus, h is defined on all positive integers and corresponds 


to the value of g(n) along the path of steepest ascent. 


Now, when A; = A; bor al? i, K Contains a function 
MiMmeomwiiecielt 15Seimportant to find the coordinate-wise 
Paemeeresteepest ascent. Without justifying its existence 
Temes pOlmlt, tie particular g(n) will be defined and a 


property of its associated function h will be noted. 


Zz 





ineomem™35-5: bet g(n) = 25 920"2) where g5(0Q) = 0 
~ n-4 

andeseei(n;) = T.Di 2 aa POL ou pos cura t 

T; peUeandi3, > 0 for all i. If h(Q) is defined for g(n q) 


as above, then h (Q) 2g (n) for all n such that 


f= (Zn)- wm, 


Proof: Let Ay yr where 
TD : ° 
ahi = U5Get) LEIS KE 
ey 


Let h(0) = 0. Finding the & — element of the sequence 
in, | is equivalent to picking the Q largest elements 
from the set Xe and letting n; equal one plus the number 
of elements picked from eri . Then h (Q) = g (nN, ). 
Eoreany oOchemen such thaty Q- (Br)y-m, 9Cn) is 
Eoura by picking the Largest n; - 1 elements of Xi , taking 
their sum, and then summing over i. Since the sum of the 


2 largest elements in Xg must be at least as large as 


vii ie 
i- (j— largest element in ) such that 


' - 
= 


ite (Z nj-m, h (2)>9(n) for aye SUEY Glatt XR = L(n,) Saye 
and the proof is complete. 


Consider L (a) weenie Ae for ald 1. 
Lins {Ace A Dm|(Tnd+ ELAS 


ws 


Loewe 1, = I 


: and suppose for the present that ite = (0 


O 
for all i. Let the set J = Pirniodt . Then 
~ Wi-t A 
Limi: {Ar A Zap {TZ + Zhn- BEE sat, 


Piemlaest cmpnession in the second set of brackets is g (n} 


in the previous theorem. 


23 





The following theorem completes the proof that 
SJemeinlate-wise steepest déscent works for L(n) and, 


RemecestOn k(n) . 


iieeuwem S24: Let ng Benderined asS@petore and let 
ie = (ne: fe O,4,:-- } ee “ih Pe). ests eet iten 


min _ min 
7 mo)! = ncW 1 


Paget. First note that theorem 3-1 and its corollary 
showed that min 1,9) Sets Come eS Ge, doe mana 


a, = Veep. + Tagieel), . Then 
en > fa, +a, Bn} ) a, = 3en)} 


leary , st Twn) Ss cur diay. suppose 

min Hcn) o< ae Loiny . Since the minimum exists, there 
- 2 - oe A ; A 

1s ann t W such that J],(n) <min Tin) : Let {= (Zn.)-™, 

Mien, trom theorem 3-3, gia) < a i = Q (Nn, ) ; 

It follows that 


JL Coe fa,+Q.Znif fa; - Oe fa,+ dy (h+re) (a,-4in)| 


2 fa, +a. (tem)} fas~ 96 \= Ly) 2 mie Liln )e 





a 


This contradiction gives the result. 


Comelilany 3-5: ne Kin) = new KM), 
Proof: See the proof of theorem 3-1. 
Theorem 3-3 assumed that fe > 0. The following theorem 


shows why this assumption is not restrictive. 
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Theorem 3-6: If Tj $ 0 for some j, then n;* Ses: 





PEoouw) ouppose that n;* > 1. Let n* be such that 


L 


= min Ltn) . Now 


(n* ) ; 
ae De obs 
Licatye far a. 2et fo eas + BP AR] 


nh = # = fe 
Let nh (n, >° © @ gy h byes Ty, Since n, Z Tees 


fa, + a,b n,) < Gar, + ae? oe ) . “Simee Jig SO 


ad 








i . i>: T.D: , 5, D: TDi Tp; 
. < « and (I, ED, + zr a Peace = ). 
Then L (nl) — Gm). “hrs contradiction implies n,* = l, 


and the proof is complete. 
For the special case of es for all i the preceding 
theorems have proven that the following algorithm will 


provide the optimal solution. 


Algorithm 


Gm Let J = {is T.>o} Wee? aleks a) - 4, SeOp. n,* ="] ton 
all i. Otherwise, set n Ot 2b). Compute L* = (one) 
EVAL De 
A. ToD: 

and go to step (2). 


and N, enc ele le eegoe cen = Cl eseeele) 


ie IDR baoe 
(2) Let j € J be that index such that mn CArL) = ane 


for all ié J. Increase n by one. If Bs - Ni: Seton; 


es n. Otherwise go to step (3). 


A 
women < bt, set 1* = L(n) andn=n. Otherwise, 


make no changes. In either case, return to step (2). 


Zo 





aera eis etound, then Q.*, ee and K* can be computed 
as shown in Chapter II (see equations (2-6) through (2-11)). 
Earlier in this chapter it was shown that for an exhaustive 
search over G for which m = 5 and N. = lb Or (1. ‘=! dene, Dante 
number of iterations required would be 100,000. If this 
avgeorithm could be employed (that is, if A; ee One 
i =1,...,5) for this case, then the maximum number of 


iterations required would be 50, an impressive reduction. 
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TV. EReAMPEES OF THE NO-STOCKOUTS-ALLOWED MODEL 


This chapter is dedicated to some examples which, it is 
hoped, will enhance the understanding of the model and will 
aid the reader in following the use of the algorithm 


developed in Chapter III. 


A. AN EXAMPLE WITH LOWER ECHELON SYMMETRY 
A two-echelon, three-activity inventory system is 


operating with the following parameters: 


Ae. > SoU ; Aa Le, 
J = 0.l/year; Lapse =e USy eau; 
C = $10/unit; D, = Dz, = 1,000 units/year. 


Tae aim 1s €o6 find the optimal operating doctrine. The 


following sketch illustrates the system. 


Qo 
Activity 
oO 
QR, OQ, 
Activity Activity 
4 Za 
D, De 


2 





The algorithm in Chapter III terminates at step (1) 


since iv Zeer Equataon (2-10) gives 


Koti, t= $21,002 (or $1,002 in variable costs 


alone), and equation (2-5) results in 
OF F-17002 Units. 


bamwat tons (2-11), (2-12), (2-13), and (2-14) yield the 


remainder of the optimal operating policy: 


Us 1.002 years; 
see = as = i; 
OV Same gy and 


Qe" 


ZOO 4suiiats . 


Had the inventory policy been derived using the standard 
EOQ model for each activity, the following would have resulted. 


Since activities one and two are identical, 


Qn: = Oe => 63.24 units; and 


Veo le =, OS years . 


This results in a time-average variable cost for each of 


the lower echelon activities of 
Wer ey | O2yyear. 


To find the time-average cost for the upper echelon recall 
Hite ly aCe nee, cilat the upper echelon must 
Supply Ce as Q,*| every T,* units of time. The holding 


and ordering costs divided by qs for the upper echelon are: 


= AD ie re & (n,-4) (Qs + Q,*) 
Cane o> = Gin Cito) 
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pamiecwO, Qs cymequation (4-1) becomes 


K= 22. + Te (nt) Q®. (4-2) 


Treating n, momentarily as a continuous variable, it is easy 
ROwsece that eCque@tion (4-2) is convex in n,. Minimizing with 
meepect tO mM, Sives two Candidate integer solutions for the 
Dermat rone propiem, my = Il and n, = 12. Of these, 
Meet l omvesmencmtowcy, COSt value for the problem parame- 
ters. Then the total variable time-average cost in this 


second analysis of the inventory system is 


Ac Dy een 4 ) Q.” + oy Ky 


n*Q 





& 


=a go Peas. 40. 


Comparing the second analysis to the first, the second 
policy yields a cost which is 41% higher than that given by 
the optimal policy for the model of Chapter II. This seems 
only natural since the model of Chapter II allows for 
cooperation among the various activities of the inventory 


system. 


B. AN EXAMPLE WITH LOWER ECHELON ASYMMETRY 

The following example has some differences in the 
parameters of the lower echelon activities and the resulting 
cost difference between models is more dramatic than that 
rience mpae VlOUSs seXanpie. 

Cons tec neameve-ceneton, three-activity inventory system, 


as in the previous example, with the following parameters: 
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A 


5 = $500; Ay= $100; Ag= $1; 


ty 


C07 anit, — £00 units/year ; 


Oi. Wegear ; pel >. = 0.05/year; 


anda, = 1000 units/ year. 


Again, since I, = I, ¢ 0; ny* = n,* = 1. In this case 


the total variable time-average cost is 
K* = $813.08. 


The optimal operating policy is: 


Ore eo an tS, 
Ory e= tes. 5 UNIS ; 
a ee = 1.478 years; and 


‘Oh O20. 15 -una ts. 
O 


As in the first example, a policy can be derived assuming 
independently operating EOQ activities and the variable costs 
can be compared. 

The EOQ operating policies for the lower echelon activ- 
ities are: 

Q,* 


Oj = 20enco muni ts ; 


ZOO utile ts. 


T,* = 2 years; and 


ie .06324 years. 
This yields time-average costs for the two activities of 


K,* = $100/year , and 


$31.62/year. 


0 


There are two alternative policies for the upper echelon to 
employ. One is to order for both lower echelon activities 
on the same cycle. This leads to the constraint 

nyT,* = n,T,*. The smallest values of n, and nz which solve 
Mis cQuat?on are nN, = foo L idee = o0,000 and ete, 
result in a policy with ridiculously high costs for the 
upper echelon. The other policy is to order for the two 
activities on independent cycles, and it is this second 
policy that is derived. The variable cost for the upper 


echelon as a function of time is 
Z(t)= Ao Intlag|+ATntlat] . 3 
o+Mtlate}+Aetntinit*] + holding costs (t). (4-3) 


Dividing equation (4-5) by t gives 


ATatlr)  A-Tat brn] 
<< 








K (t)= i + Qve. helding costs (t) 2 (4-4) 
As t grows large, Int ates! ager a large integer; and 
—— 
Int lee] 2a) &a | 4 
1 t f. () (w.T, : “Th : 
Similarly, 


Tat Le eek 


nT” 
Average holding costs, by an argument similar to the one 


above, stabilize at 


ie C (n,- 1) Qs FA ine (n2-1) an 
a a 


Then lim K(t)= Ke where 
+00 


oun 





ef , Ae ae 1) Qr T.C(n,-1) Q> 


nie ee oe a eS) 








= K(n, >My). 


Equation (4-5) is convex in n, and Ny if n, and ny, are 
treated as continuous variables. Minimizing K with respect 
to ny and nN, gives n, = 1.58 and n, = 15.8. Four alternative 
solutions for integer values of (n,n) are (1,15) 2 (2 eis) 
wiero). and (2,16). Of these four, the choice that offers 
the smallest value of K is n,= 2 and n,g= 16. For the 


mpper echelon alone, then, 
K(2,16) = $1,193.45. 


Adding to this the time-average cost for the lower echelon 
activities gives a system time-average variable operating 
Gest Ofesl, 325.07. 

The second operating policy yields a system time-average 
Geet whiren 1S 635% higher than that for the optimal policy 
fomeaemwe mogel presented in Chapter [1]. With asymmetry in 
Mire problem, allowing the activities to cooperate pays even 


greater dividends. 


C. AN EXAMPLE WHICH EMPLOYS THE ALGORITHM IN CHAPTER III 
This final example shows more details of how the 
alcoriénm of Chapter III works. 
Consider a two-echelon, three-activity inventory system, 


as before, with the following parameters: 


SM 


A 
O 


I 
O 


C= s10/unit; Dy 


5s A, = A, = $5; 


ay ercat ly =) 0 l/yedr; 


100 units/year; and 


De = 1000 units/year. 
Domne ene), uin) = 316502)" Then 


1650-50 
25 


Nee ime pooh 00 sap +1 = 6. 


ae ce Pei = Oo mand 


eamee for 4 = 1, 2 


I, Dy =-5/2 @and 
n(n,+4) 


I, Da. 
Naln.+4) 


set ng = 2. Now n = (1,2) and L({n) = $1700. 

Then n remains (1,1) and L* remains $1650. 

The table below shows the progress of the algorithm. 

In this example ny4* = ng* = 1. The optimal policy could 


then be computed using the equations at the end of Chapter Il. 





it00 | — (1,4) 






1650 Ci) 


-% 2 
ie 
“EQ 


SCSI Cen = 
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V. THE BACKORDERS-ALLOWED MODEL 


A. MODEL FORMULATION 

This chapter extends the model developed in Chapter I 
by allowing the lower echelons to run out of stock and to 
backorder the shortage until the next order is received. 
For each unit backordered a penalty cost is assumed to be 
suffered. 

The sketch below shows the net inventory during one 
cycle for one lower echelon activity. For activities one 


through m a new variable, s the backordered quantity for 


1? 
one cycle, and a new parameter ue the shortage penalty 


cost, are introduced. 


Wo 1) 





The total costs for one cycle are given by (5-2) 


Fz Ag+ CUnQ+ TCIR EE Qe Dah 


+ © (holding cost for lower echelon activity i) 


+ © (shortage penalty cost for lower echelon 
leteniny aa tye i ya Cs) 


a5 


The holding cost is the product of I.C, the unit holding 
cost per unit time, and the area under the on-hand inventory 


curve over n, pemgaeds. That 15: 


s.)" 
ne (he) £(@r-8)T= nT e GE 


The shortage penalty cost is the product of Ts ess timeee 


penalty cost, and the total number of backorders incurred 


over n; periods. That is, n;1;S;. Thus, equation (5-2) 
becomes: 
Cnr- LIN: 
Z Z A. r yn Ai eG bn Qi free)... ie Q. 
ny (Q.-S.)* } 
SL, IC ie es Le (5-3) 


Dividing equation (5-3) by i gives the time-average 


cost for operating the inventory system. 





eee Cm DiQy ni <i 
K = aon 7 Oe et sees, and, +Z MAL 


ne De Qi s,* Wn D| Ga ND, 
= a —_ ° ee + ne —_— = 
De 2n dD, oa + 2n,DiQ ds n,Q, ® és 4) 


Equation (5-4) will be denoted by K(Qq 55, 5+-+-5S5,.My,--- yn). 


Peo mien OPERATING @OCDRINESFOR GIVEN n ,...,n 


a 
A full statement of the problem faced in minimizing 
time-average system cost would be: 
minimize KG: Qn... Tray Si, «= Sia, ---) Mim) 
subject to Qt LN@Qi and Qiz= Dil for it 44-53 
n.2 4 awd integer feat eee 
G20 for istynm 
erssO mcr > 
and nT =e nyT, = To For te dy ym and jodie. 
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Equation (5-4) presents K with the equality constraints 
Gibecaay sdtaciices  fuli justification for this simplifica- 
tion is not presented here, but the reader is invited to 
consider what information the Kuhn-Tucker conditions would 
give in the case of the equality constraints and how this 
would then be used to simplify the other first order condi- 
tions. Consider only the problem posed in minimizing K 
given by equation (5-4). 

minimize K(Qq Sy a--+ oS My +--+ Nd 
subject to Q, 7 0; 


ON<aasey'< EES Oi: and 


a n, Di 
n; 2 i@odiad “1k per . 
Disregarding the variables N4,-+-5N, for the time being, 
Wetice that the Constraint set is convex. If K itself is 


convex, then the entire problem is a convex programming 
problem with a global optimum. The issue of the convexity 


of K will be discussed later. 


2 A A 
Let Q4sS45---55) denote the values of QoS eee oS, 
which minimize K for a given choice of Ngy---5T. Let 
he 2 Se = »S*.1,%; » Ae ale be the optimal values for all 


the variables. By arguments similar to those given in the 


proof of theorem 2-1 the following can be shown to be true: 


e a # 
inte Come moe COges mer Sattar eyrin = (CQL Si, ...5Sem 5 ttynsy Mme) 


e 


A A 
Now eee can be found by the calculus. 


(Cs, oS pa oe) is a global minimum for given Ny s--+-5N, 
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Picante amaldieaGderivatives of K with respect to 

Q4 Sq o-+ 05S, evaluated at OMB 5 Sy isos oN, are zero 

and the Hessian of K is positive definite (so K is convex) 
OG wilews peer. ted Nyy++-,n, and for all feasible values of 
Q42S4 ee 29S) There are several cases, however, in which 
this approach to solving the problem fails because either 
Picmpontdal derivatives Cannot be set equal to zero simul- 
taneously for some aS 0 and some Sr such that 

Os A $ Sp G @, for each i, or the Hessian of K is not positive 
definite. 


Taking the first partial derivatives of equation (5-4) 


Vereie respect torQ, and $5: 











ao K . a AoD, Gr Din: | nv AD, 
TQ, ~ n,Q.” + [1.0% an D, ~ Q? 2 n, (S25) 
T.c n,D:z aD | a € s° na) | TUS.N.D, 
u L AnD, 7 qe Lo abe ~ Q2 De n, peel 
PK Toe se Tee by 
55,° LO + “pe 7 mea (5-6) 


setting equation (5-6) equal to zero and solving for S; 


Maielas : 
D: Gin. C ~ TiniD,) 
i ear a : (5-7) 


Substituting this result into equation (5-5) and setting 
(5-5) equal to zero gives a quadratic equation for Q 


in terms of the system parameters. 


B, ae mm en = O where (5-8) 
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n, D. I. Cc 


B =m ae anv Di - 4 
74 Cni-')n, Di mT.C D: and 
Got, ee vn oti 
T an, Dd, - 2 ZN Dy ’ 
AeP, n- Di Tt. D, 
i QAn,T.C 
- ee 
2 Cng-t) My, Di ibs nd & 
I.¢ Sm 2n.D, . AnD, : 


Depending on the values of By, and By, the method for 
finding Qy Sq oe eS, varies so that a case by case analysis 
will have to be made. Note that it is possible that B, can 
beeTllegative, Zeno, Or positive, depending on the sign of 


2S muNnerator, 


Peed, nm. AiD, n: De Ti D, 
ps L " 7 An, Lic * 


" (5-9) 


Note also that if n; = 1 for all i then By= 0; otherwise, 


Bia< O. Six cases will be considered below: 


CASE By Ba 
ONE = 0 = 0 
TWO <n) > 0 
THREE = 0 > 0 
FOUR < 0 = 0 
FIVE < 0 <0) 
SIX = 0 < a5) 


In the analysis of the cases, the following property 
of K will be valuable. For a fixed, positive Q,, K is 
positive definite everywhere. for Sqye-eaS)- ‘1 . If some 
of the oe = Poni emo eCICT ro. pOoteave dctinite eveny — 
where in the remaining S$; 7 0. To see that this 1s true 


a 
consider 35,98; 


oo 


aK ioe nD, 
oS 7 nN,.Di Qi for i = us and 





zs O oie elle. |). 
Hence, the Hessian of K in any non-void subset of the S;'s 
is a diagonal square matrix with all of the diagonal terms 
positive. This implies its positive definiteness and the 
resulting convexity in the S,'S. Additionally, when the 
value of Ss that renders 2= = 0 for each i is substituted 
peo Kegiven by equation (5-4), Kk is then convex in Q,. 

This method for finding the minimizing values of 


Qqy,Sq,---, ands, ise sttiteds im Seck1on D. 


Cree ONE: "Ey = 0 and By = 0. 

In this case equation (5-8) cannot be used to solve 
for Q,4 since it yields no information. Let Q,> 0 be 
arbitrarily large but fixed. Then S; 1s given by equation 
(5-7), and the positive definiteness of K in Sqyeee5S) 
implies that Messe provide a global minimum for the 


given Qy. Now consider K as a function of Q, at the value 


A A * 
Sqo+++5S)- Its value 1s 
= GoebD. + 2am. Di . 
A 
pedecs) ots Iidepem@ent of Q,, Q, could be any positive value. 
Once it is selected, S could be computed using equation 
(5-7). 


A 
A problem arises, however, if Q4 is chosen such that 
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Qs S teigm ae ae 


since equation (5-7) would give S5 < 0 for some i. The 


problem is resolved by considering the following two state- 


Mewes sor tne problem. 


Aw 
Minimize K(Q,,S,,.---5S,) ; 


- max 1.D 
Q, 


given 4cvim Toe 


A 
Minimize K(Q4 5S, peer aS) je 


max Tid. 
given Oh <qeiem Te ; and 


S; > 0 for all 1. 


Since the second statement is a restriction of the first, 
the best answer to the second can be no better than the best 
answer to the first. This implies that ee srs Oss 
Oy < 4siem a is greater than or equal to CLIi+ ZU D: 
and that there is no advantage to choosing 0, that small. 


A A pra ee Dt 
Thus, Q, can be picked as long as Q, 2 qicm Tec: % 





A 
and S5 is then computed using equation (5-7). 


CASE TWO: By < 0 and Ba > 0. 


A 
iments «Case sO Cdn be computed using equation (5-8); 
A Bz |% max Mi Cids 
Mame Liya0 04 = |- an a bOvedcds tha t Q,2 isigm nN, Tic 9 


A 
$5 will be laa ca, Sakae Pome dll Vv. ana@=etne problem willbe 


max MED, 


solved. If 04 < Agew nT. C , then oF <0 forssome 1 


max "mid 


MOUs result.. Thus. for Fall n tic 


, this method 


Ceesolielen 1s) Invartiad. 
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nv id, 


~ Ma x 
Suppose that Qy 7 4cigm N oa oon and fixed. Solving 
tO I; Br Se that minimize K results in 
\ ox n,Cm.-1) D; 
Ks a. | expression 5-9)| +CZdi + 7, Tid: + Q, {roc yc Bal 


max Mm TDi 
Additionally, suppose that 0, <iscgem nic « 


Then K (Ors SCA, SQ) < minimum K subject to 
Qi = O14. In fact, as Q, decreases from Qy towards on, 
K decreases. Let 
ma x n.TD, 
and let Q, decrease from 0, Pome ue e,- Minimum K is taken 


on where S5 = 0 for at least one value of i. Let 
Yy =} i Ss, = 0 { : and 


. mts, > 0} 


Now as Q,; is decreased, the minimization of K given Q, < e 


Y 


al 
becomes 

minimize kK 

Si: cz Y2 


subject to Q,4 fixed; and 


S.9= U0 ror 1:2 Yq . 
ai 


ive fUmetronal form of minimum K is: 


a Ae D es nip, T Di 
Ke a (SP + 2 Bee Roane | + €EP 


Q an 
.-ln,D T.c nD: 
+ cars Di sis @, {Tcl ae D, tL ID, 2nd, \ (5-10) 
ae 


The first bracketed expression remains positive and increases 
as the number of elements in Yy decreases. The second 
bracketed expression also remains positive and increases as 


the number of elements in Y, increases. 
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Equation (5-10) remains a valid expression for minimum 
K in SqaeeeaS) given Q, fixed until Q, is decreased to the 
point where equation (5-7) implies that another Ss; = 0 i OG 
some i. Then the functional form of minimum K changes, 
because the sets Y, and Yg change. Once Y, and Yg have 
DecheM@citled tO GeEStrice, any of the S,'s from turning nega- 
tive, equation (5-10) is again valid. This method of 
recursively formulating K remains valid until Q, has been 
decreased to that value which requires Y9g = @ and 
1S fi,...,mf. Then the problem reduces to the one given 


iioacaapteredt , S$; = 0 ora Liga. 





As 1 ranges from one to m, let Cy > Cn> 2+ VED be an 
_ a Ni. M.D, ‘ 
ordering of the ratios Seren 0) that 15,.-lhet 
Wee G 

Htc 8 % 

ey = max ee Ue 4,..um} ) 
Nid, ea | 

= eae . = eee ; Qna 

€, = max [ae Came = tyrone ( 
Ni. D: a m| 

e. = max CGem ) C= be-°5 : 

P ec ut ; 


am | 
If Qy < e4, the optimal values of Qs»Syo+--, ands, are those 
which correspond to the least value of K from among the 


sequence of problems stated in the above two paragraphs. 
The first problem in the sequence is to 
Minimize kK 
SD) 
Subject to Siz O for ve Ns 7 and 
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where Yq and Yq, are the first Y, and Yg formed. 
For this oe 


+ E, + Es. Qs where 


Eee > as y, Midi 3 "and 





LLY2 
.-t)n, De ae nD. 
E 2 = ras ZS aE D, + ny oD 4 a 


Now K is convex in Q,, and, hence, minimum K corresponds to 


ae: ta 
Ones {== if ae { €:| £$<€, 3; otherwise 
0. 7 e. Soh that K Cee) = a Ke. ° 


thse pneccaune 15 tollowed p-1l times for the p-] inter- 
vals for Q,, and for each iteration a value for minimum K in 
Og is generated. The value for K where oa QO for all i need 
not be considered since the last of the sequence of p-1l 
problems examines K(e,)- From among these p-l alternatives, 
pick the least value for minimum K. This will give a feasible 
value for on from which feasible Ba can be computed for all 
iepe Obethat iteratden. Of course’, oe Umer al) ker 


Ot eghat 1heratrons 
CASE THREE: Ea Oe act B, > OL. 


In this case equation (5-8) cannot be used to solve for 
A 
Cho dee -& | PiOeeacuiicdnme Fete. > 0 be aGbmeman ily 


A e ° . 
large but fixed. Then S, is given by equation (5-7), and the 
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positive definiteness of K in Sqyeee 5S) ip LTes hat aon 
ipieced. Q.., eee Sk provide a global minimun. 
Now consider K as a function of Q, at the point 


(eS ye 


m 
ee = [expression (5-9)) + CLD + LUD , Cs- in) 
\ 
For Qy 2 e€,4, equation (5-11) is valid and minimum 
f= liramkK = C TDu+ Ze Di . 
Q,—> 00 
But suppose that Q,y < e . Then equation (5-11) is no 


longer valid, and an analysis similar to that performed in 


CASE TWO is required. The only modification is that Ea 





Leo: ae . 
reduces to E3 = Des 2 = The minimum K is picked from 
Y, Dy 
among the p-l alternatives. If that value of K is less 


A 
than CLD.+ Zid, then Q, oe ares se are the corresponding 
values for which that K applies. Otherwise, Qy =®© and 


A 


Ss © for all i. This suggests it is best not to go into 


business. 
CASE FOUR: Bw< 0 Sand Gba.- 0; 


In this case equation (5-8) cannot be used to compute 
04 Since it implies that O4 = 0 which in turn would imply 
S. < 0 for all i. As before, minimum K for a given large Q 
can be found when sm is determined by equation (5-7). 
then K = C2 Dee+ 2D, + Q {Toc 1 OAD: | . Since K is 
Seley einechrcasimeo= in Or, Kk 1S minimized as Q 4 1S minimized. 


But this argument is only valid for Qy 2 e,. Hence one 


candidate for minimum K is 
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Chea ine. 
K= CLD + CEDi+ & [Tot Land. 


(CorenleZ)) 
For Qi < & > 5; 8 0 for some i; and the formula for minimum 
i “enanges . 
Since Ss. = QO for some 1, B increases; and, hence, 


Be > 0. Now CASE TWO applies. The minimum K under CASE 
TWO conditions is computed and compared to equation (5-12). 
Peierine (soe) 15 less than the CASE TWO solution, Q4 = e, 
and = is computed with equation (5-7). Otherwise, the 


CASE TWO solution applies. 
aoe FIVE; By, < 0 and B, < 0. 


In this case equation (5-8) cannot be used to compute 
: B. 1% ; , : : 
Ono 1mee = a 1s imaginary. As before, for a given 


large value of Q4 the minimum value of K is given by: 

k= & [expression S-a)] + CID. +L Edis Q {rez SOY 

(5 -sles) 

Since expression (5-9) is negative, K is minimized for min- 
mat Oy. Bat equation (5-13) 1s only valid for Qy2 Ee. 
Hence, a candidate for minimum K is the K given by equation 
(5-13) evaluated at Q, = e,. 

For Qi < €, construct Y, and Y, as was done in CASE TWO. 
Here, K is given by 
nope ao: 


A) 
C= a ieee - L Ontee | +c DD: 


(ni-L)n,Di Tec we 
2 ee? + Q, ire 2nd, - u An, Di 
Y2 (5-14) 
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Now, if oe ti, 
: 7 ne A.D Now ey t 
Be ep ee Es 2c | 


| 7A (5-15) 


remains negative, then the next candidate for minimum K is 
K given by expression (5-14) evaluated at e,. This analysis 
continues (evaluation at e:, would be next) as long as 
epxression (5-15) remains non-positive. If expression (5-15) 
turns positive at some iteration, then the procedures of 
CASE TWO apply for the remainder of the analysis. Then 
minimum K is picked from among the p alternatives, and 
Q..8,.---, and 5, are the values that pertain to the K thus 
paeked . 

Although this search is sufficient to find minimum K 
it can be streamlined somewhat so that fewer alternatives 
need be considered for minimum K. As K is transformed from 
iteration to iteration and indices migrate from Y, to Y,, 
iSisecontanuous in Q,-. That means that K for two successive 
iterations when evaluated at the common boundary point is 
the same. While expression (5-9) remains negative we have 
seen that minimum K in Q, is taken on for Q4 equal to the 
left hand or smaller boundary value. Thus, as K changes 
from iteration to iteration and as long as expression (5-9) 
remains negative, K continues to decrease with decreasing 
Q4. Accordingly, past values of K may be discarded as 
potential minimum K since they are greater than the present 
Ke 

The first value of K that should be retained is the last 
value of K for which expression (5-9) is negative which 
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equals K evaluated at the right hand or upper boundary value 
BOENOv LOR Wiich tleratron expression, (3-9) 1s positive for 


the first time. 
CASE SIX: B, = 0 and B, < 0Q. 


This case is the same as CASE FIVE except that the term 





Via Habe : : : 3 
Q, [re 2 SOR | is missing. If expression (5-15) 
turns positive, then this problem falls into CASE THREE 
SMeeic tat Eee 2s ELC Di The appropriate 
3 Y, AD, 


A A 
values of OWN >) encase are found as they were before. 


C. THE SEARCH FOR OPTIMAL Ny,---,N, 

Finding optimal Ny>---on, involves many calculations. 
The first problem that must be addressed is, as in the 
no-stockouts-allowed case, showing that only finitely many 
peincs need=to be Considered. In general this cannot be 
done since, without some restriction on the parameters, as 
the values of N4,++-,n, range, the solutions for 
Q4,$,5---, and $ found in Section B cannot be controlled 
and the value of minimum K cannot be predicted. However, 
if the assumption that 2AUiC Se ebe fepe GE 1 1s 
made; the hyperbox G can be constructed. 

Consider the following two problems: 

(1) minimize K(Q,,5,;... 9S Ng ress yn) 


subject to n; & Il integer; 


Q42 0; and 


n,Di fa> 3 
ae for allt 3 and 





oS 
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(2) minimize K(Qt S54 o+++5S5 My s--- oN) 
subject to n, 2 l integer; and 
Q, > OF and 


iD, 
Ss ape fOr yavlelemiee 


The first problem is the minimization problem that CASES ONE 
through SIX address with Ny,-+-,N, fixed. The second problem, 
eivem that 2A:1-C 2 ee has as a solution S and 0, 
given by equations (5-7) and (5-8) for a given Nyy---,N. 
The second problem is also a relaxation of the first problem. 
Let Ky *(my,---.n,) be the optimal solution to the first 


problem for a given Ny y+++ nN Let K,*(my,---sn,) be the 


oe 
corresponding solution for the second problem. Then 
* + 
K> Ga, n ) S Ka (nia eety Mmm) 
ite toliliowing tneorem Shows that n* = (1,...,1). 


Theorem 5-2: 2A;1;C 2%, D, for all i, 


eicieneees (lo. , 1). 


Proof: The following argument shows that, if 
for all i, then whenever n, 2 2, regardless of the choice 


‘ope ns» 


KF << ae Cle a Kg* (ny ,--- MJ) $ Ky eee 


= x - x). * = 
where K* Ky (ny peee oT, )3; and, hence, n; ll fon eaehwe: 


m 


Bloc eHEUting Q4 given by (cquation (5-3) and S; given by 
equation (5-7) into equation (5-4) for K yields equation 
(5-16), the solution to the second problem for given 


es 
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VV 
Ka “(n= By+ [Been Betar4 Bey (Bet B+ BO)] | * 5-16) 
where 


Bons ?_ CD; +), uD: ‘ 


_ Cn.-') Di 
Z n, (2MTEC -T2D:) 
esa) = *o*+ 2 Gn 
wee 
By ma AG 9 


es (Dn Ai)* ; and 


a min, (1:C - 271 D,) \% 
oe” eae) Ce 


Note that if n. 2 oO SOM et (Onn rc U5 eee CONdhet On 
that makes this solution for K valid), then B3, Byg(n), Bs(n), 
Bg, Bon )e@manesBe(n) are all positive. Then, for K given 
by equation (5-16), 
kK.) > Ba. 
Consider now K(1,...,1). Since B, > 0 and By = 0, 


Gist THREE applies. Recall that for CASE THREE 


K(Q.)= B, 


min {(Q.) s lim min 
Siy95m Qe Suse 


It follows that whenever Nn; 2 2 for some i 
* * e 
Gale) evap Catone (ny mrcemiggicnen) 
Thais Wompletes the pwoof. 


) 


Throughout this chapter the various efforts that have 


Deeps CONVEXITY OF K(Q4q.S5+--55, 


been made to solve for minimum K have disregarded the notion 


Pictekenay not be convex an Qy,si,.-.., and s_. In general, 


20 





ier 1S ono t. 


One way to show convexity is to show that the 


Hessian of K, 


mesttivwe cderinite. 


Now, 


WU 
~ 


XK 





Vv 
£) 


Y, 
x 


“Pp 


| 


VV 
~Y 


35,22, 


Let Ch; 


LOL sis. 5 , 


28, out of 100, non-zero entries). 


(hy 


case for m 


WU 
“A 
KX WV 


) 
(h, 


;? for a problem of arbitrary size is difficult. 


the matrix of second partial derivatives, is 


The Hessian of K is 








°K O*kK Ork 
90,95, 30.252. a ts 3Q,35m 
sk otk 
een ostos 8 | i 
‘ 
ee ee ee re 
aS 
ZAod, | paid i . 7 SEMBES , 7 HSM 
n,Q? a 2a. n, ) 
-T,Csind tind . 
Qin Di AQ 
Ae SOs D, sean ee 
Vig DD. Q, For oe ane 
O Sor Lts. 
, ot I< 
be the Hessian of K. Since S55 = © 
Oxy 
>? is sparse (if m = 9, then (hs) has only 


However, the analysis of 


The 


J] illustrates this. 


ol 





suppose m= 1. Then rae has only four entries and is 
posrtilve definite provided that h,,? 0 and as haa hin {> 0. 
For Q,> 90, hy > 0. 2 1 sos : cis? tia | 
ICA ° § as es 
hihas- ha rs By {Ros +A, L,C+ +I¢r o Sn ue . 


“Dt Dat 
Peeae oo, tas dees to 


2 {Lcd , are-n}| 
Q,* a 


which is negative for M1, sufficiently large. Hence, (hj ;) 
is not positive definite for all feasible values of Q; and 
Sy; and K is not convex. 

Although the convexity of K is a sufficient condition 
to claim that a feasible stationary point of K is a global 
minimum, it 1s not a necessary condition. The analysis in 
"ection B showed that if K was first minimized over 
Syy+--, ands, (K is convex in S4,--., ands_) and if the 
resulting cost function was then minimized over Q, 

HAY 5) p<: - 8) 1s convex in Q4), a global minimum could 


be found. The following theorem formalizes that approach. 


Theorem 5-3: Let 


CQ.) = KCQSiyendm) = pimimpe KCQuSty 95m), 


§ mr 


and let , 
K(O, ) - min K(Q,). 


Qi 
A A 7 ~ 
Then K(Q4*,5 *,-.-55,%) = K(Qi)3 Qa* = Q5 and s;*=s;(Q) 


homme acn 1 . 
Proof: Since K(Q*,5¢*,---55,%) is the global minimun, 


K(Qy*,s *,.--550*) § K(Q). 
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Now, 


K (Q4) $ K(Qq .S4*5-++55 


and, hence, 


A A 
K(Q4) 


a ae 


m 


1A 


K(Q,4*,5, *,. Be a) : 
Therefore, 
m= K(Qi* see eee) 

Bemelily, 1t follows that Q,* = Ox ; ands: “ = S(O) oo 
each i. This completes the proof. 

This method of finding the optimal solution is used 
throughout CASES ONE through SIX. It suggests a weaker 
sufficiency condition for (Q4*,54*,.-.,5,%) than convexity 


of K. This condition is proposed in the following corollary. 


A 
Gomollary 5-4: Let S; be the value of S; which minimizes 
aA 
K(Q4.S45+-+25,) for a given value of Q. Then K(Q,) = 


K(Q4 5S, (Q4) +--+ 5, (Q)) 1s convex. 


E. CLOSING REMARKS 

Uhasschapter has offered the development of a backorders- 
allowed extension of the model presented in Chapter II. 
Section B derived a method that can be used to find candi- 
date points for a local minimum to K. Section C showed that 
n,* = 1 provided that 2A-15C Zee D. for each 1. Section D 
showed the solution technique discussed in Section B indeed 
provided the global minimum. In all, the general backorders 


problem has been solved except for the problem of bounding 


2 
phewsearch in n when 2A.1.C< 1, D; for some i. 


oo 





VI. CONCLUSION 


A. SUMMARY 

This paper has presented two deterministic multi-echelon 
inventory models. The assumptions of the first model led to 
the development of an extension of the EOQ (no-stockouts- 
allowed) model. Through the presentation of a handful of 
theorems, it was shown how the optimal policy could be found 
either by exhaustive search in the completely general case 
or by a rapid steepest descent search in the case where the 
ordering cost is the same for each lower echelon activity. 

From the author's point of view the three most inter- 
esting aspects of this paper are the assumptions for the 
Steepest descent algorithm, the effect of i = [. = t on 
SorimMal mM im the mo-stockouts-allowed model. and thesgenerall 
method for finding optimal K(n) in the backorders model. 

The first point of interest is the way in which the 
assumption A; ="o;, fOr Cachei alteetsmuncmsiea) chtOm opr iiiam 
n in the no-stockouts-allowed model. That this assumption 
streamlined the search is of some interest mathematically, 
but it is of additional interest from an economic point of 
view. In military applications of inventory theory it is 
usually assumed that the cost of ordering or requisitioning 
is constant for any one class of activities. For example, 
among Navy stock points inventory "pull" item orders have 


roughly the same cost. In the case of afloat activities 
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it was long assumed that the cost of submitting and pro- 
cessing any requisition was $7. The exact figures are not 
important here, but the principle lends additional justifi- 
cation for the assumption. 

The second point of interest is that the observation 
that It = I.-1,¢0 implies that n,* ee LS Say Serie eet 
cost to the system of holding goods at the upper echelon is 
weeleadst as great as the cost of holding them” at that par- 
ticular lower echelon. In this case, the solution shows 
mace the Optimizaeion Strives t@ minimaze the holdings of 
ehe upper echelon because of the relatively high “opportunity 
Gost’ presented. This notion seems to support military 
decisions to provide central procurement activities with 
an unencumbered "open to buy," because the "opportunity 
costs" for alternative decisions are high. 

icmeiiudmpoint Of Interestaisesemctly a mathematauca 
one. Even though the cost function in the backorders-allowed 
model is not convex in all the continuous variables simul- 
taneously, theorem 5-3 and its corollary show a global 
optimum can be found by decomposing the cost function 
minimization problem into a sequence of two convex-program- 
ming problems. Thus, for this problem, a weaker sufficiency 


condition than the standard one of convexity has been found. 


B. RECOMMENDATIONS FOR FURTHER WORK 
Specifically, there are two pOdnES) Mn tints paper sundae 


should be re-examined and, perhaps, extended. The author 
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feels that an optimal search plan, such as the algorithm in 
Chapter III, can be developed for the general no-stockouts- 
allowed model. For the backorders allowed model it was 
: 2 : k = 
shown that if 2A-1,02 T: Ds for all i that n; 1 for 
all i. A bound and search scheme should be developed for 
2 
mie case where 2A.1.C< UC. D. for some 1. 
ea! ala eal! 

A general area suggested for further work is to extend 
the backorders-allowed model for the case in which the 
backorder-penalty cost may be time-weighted. Additionally, 


a "lost sales" model could be developed. 


C. CLOSING REMARKS 

In appraising the value of this new model it should be 
remembered that it represents an attempt in basically a new 
direction. Simple results usually forerun more complex ones 
which are often more directly useful. However, just as the 
standard EOQ model serves as a good first approximation to 
Heres complex deterministic or stochastae sangle activity 
models; so perhaps this model may be found similarly useful 


for multi-echelon models. 
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